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O , Abstract 

(N 

We study the Hamiltonian dynamics and spectral theory of spin-oscillators. Because of their rich 
j^ I structure, spin-oscillators display fairly general properties of integrable systems with two degrees of free- 

dom. Spin-oscillators have infinitely many transversally elliptic singularities, exactly one elliptic-elliptic 
singularity and one focus-focus singularity. The most interesting dynamical features of integrable sys- 
tems, and in particular of spin-oscillators, are encoded in their singularities. In the first part of the paper 
we study the symplectic dynamics around the focus-focus singularity. In the second part of the paper 
r tn . we quantize the coupled spin-oscillators systems and study their spectral theory. The paper combines 

j-y^ I techniques from semiclassical analysis with differential geometric methods. 

t^ ! 1 Introduction 

Coupled spin-oscillators are 4-dimensional integrable Hamiltonian systems with two degrees of freedom 
constructed by "coupling" the classical spin on the 2-sphere S^ (see Figure [3TT] ) with the classical harmonic 
^ ' oscillator on the Euclidean plane M^. Coupled spin-oscillators are one of the most fundamental examples 

^ ■ of integrable systems; their dynamical behavior is rich and represents some fairly general properties of low 

■^ i dimensional integrable systems. The goal of this paper is to study coupled spin-oscillators from the point of 

O ' view of classical and quantum mechanics, using methods from classical and semiclassical analysis. 

l/^ ■ A 4-dimensional integrable system with two degrees of freedom consists of a connected symplectic 

O , 4-manifold equipped with two almost everywhere linearly independent smooth functions which Poisson 

commute, i.e. two smooth functions on the manifold such that one of them is invariant along the flow 
of the Hamiltonian vector field generated by of the other. The most interesting geometric and dynamical 
features of integrable systems are encoded in their singularities, i.e the points where Hamiltonian vector 
fields generated by the functions are linearly dependent. Around the regular points, the dynamics is simple, 
C^ ■ and described by the Arnold-Liouville-Mineur action-angle theorem. As we will see, the dynamics near the 

singularities is in general much more complicated and depends heavily on the type of singularity. 

Let us explain the construction of coupled spin-oscillators more precisely. Let S^ be the unit sphere in 
M^ with coordinates (x, y, z), and let M^ be equipped with coordinates (n, v). Let X, p > be positive 
constants. Let M be the product manifold S"^ x M? equipped with the product symplectic structure XL0g2 © 
pojQ. Let J, H : M — )• M be the smooth maps defined by J := p{u'^ + v'^)/2 + Xz and H := ^ {ux + vy). A 
coupled spin-oscillator is a 4-dimensional integrable system of the form (M, Aa;52 © pujQ, (J, H)), where 
ujg2 is the standard symplectic form on the sphere and ujq is the standard symplectic form on M^. 

The singularities of coupled spin-oscillators are non-degenerate and of elliptic-elliptic, transversally- 
elliptic (both of these types are usually refeiTcd to as "elliptic singularities") or focus-focus type. They have 
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infinitely many transversally-elliptic singularities (along a piecewise smooth curve, as we shall see), one 
elliptic-elliptic singularity at (0,0,-1,0,0) and one singularity of focus-focus type at (0,0,1,0,0). The 
J component of this system is the Hamiltonian (or momentum map) of the S^-action that simultaneously 
rotates about the vertical axes of the 2-sphere, and about the origin of M?. The H component is given as 
follows. Using the natural embedding of S"^ in M^, let vr^ be the orthogonal projection from S^ onto M? 
viewed as the z = hyperplane. Let (x, y, z) G S'^ and (n, v) € M^. Under the flow of J the points 
(x, y, z) and (n, v) are moving along the flows of z and {v? + f ^)/2, respectively, with the same angular 
velocity. Hence the inner product {iTzix, y, z), {u, v)) = ux + vy = 2H is constant and commutes with 
J. 

Because H does not come from an S*^ -action, coupled spin-oscillators are not toric integrable systems 
- they are what now is called semitoric integrable systems, or simply semitoric systems. Semitoric systems 
form a rich class of integrable systems, commonly found in simple physical models. For simplicity, through- 
out this paper we assume the rescaling \ = p = 1. The statements and proofs extend immediately to the 
case of A, yO > 0, but we feel that the notation is already sufficiently heavy so we shall avoid carrying these 
parameters. 

Semitoric integrable systems 

Our interest in semitoric integrable systems was motivated by the remarkable convexity results for Hamilto- 
nian torus actions by Atiyah HI, Guillemin-Sternberg |[T4l . and Delzant ||6l. Despite important contributions 
by Arnold, Duistermaat {81, Eliasson lITOl . Vu Ngoc Il23ll25l . Zung ||29l and many others, the singularity 
theory of integrable systems from the point of view of symplectic geometry is far from being completely 
understood. As a matter of fact, very few integrable systems ai^e understood. The singularities of these 
systems encode a vast amount of information about the symplectic dynamics and geometry of the system, 
much of which is not computable with the current methods. 

This singularity theory is interesting not only from the point of view of semiclassical analysis and sym- 
plectic geometry, but it also shai^es many common features with the study of singularities in the context of 
symplectic topology ll20l[T6l . algebraic geometry and mirror symmetry (see |[T3l and the references therein). 

The coupled spin-oscillator is perhaps the simplest non-compact example of an integrable system of 
semitoric type. Precisely, a semitoric integrable system on M is an integrable system J, H ^ C°° (A/, M) 
for which the component J is a proper momentum map for a Hamiltonian circle action on M and the map 
F := (J, ff) : M — )• M^ has only non-degenerate singularities in the sense of Williamson ETI . without 
real-hyperbolic blocks. This means that in addition to the well-known elliptic singularities of toric systems, 
semitoric systems may hast focus-focus singularities. 

Semitoric integrable systems on 4-manifolds have been symplectically classified by the authors in 
|[T8l[T9l in terms a collection of five invariants. While conceptually they are more easily describable, some 
of these invariants are involved to compute explicitly for a particular integrable system. The most diffi- 
cult invariant to compute is the so called Taylor series invariant, which classifies a neighborhood of the 
focus-focus singular fiber of F. This invariant, which was introduced in ll23l . encodes a large amount of 
information about the local and semiglobal behavior of the system. Focus-focus singular fibers are singular 
fibers that contain some fixed point m (i.e. rank(dF) = 0) which is oi focus-focus type, meaning that 
there are symplectic coordinates locally near m in which m = (0, 0, 0, 0), a; = d^ A dx + dr/ A dy and 
F = F{m) + (x^ + yri, xt] - y^) + 0{{x, C, V, ??)^)- 



Dynamics and singularities of coupled spin-oscillators 

The coupled spin-oscillator system has non-degenerate singularities of elliptic-elliptic, transversally-elliptic 
and focus-focus type. It has exactly one singularity of focus-focus type. Near the focus-focus singularity, 
the behavior of the Hamiltonian vector fields generated by the system is not 27r-periodic, as it occurs with 
toric systems. 




Figure 1.1: Singularity of focus-focus type and vanishing cycle. Topologically a fiber containing a single 
focus-focus singularity is a pinched torus. 



Loosely speaking, one of the components of the system is indeed 27r-periodic, but the other one generates 
an arbitrary flow which turns indefinitely around the focus-focus singularity and which, as F tends to the 
critical value F{m), deviates from periodic behavior in a logarithmic fashion, up to a certain error term; this 
deviation from being logarithmic is a symplectic invariant and can be made explicit - it is in fact given by 
an infinite Taylor series (5)°° on two variables X, Y with vanishing constant term. This was proven by the 
second author in |[23l . The goal of the first part of the present paper is compute the linear approximation of 
this deviation. 

Theorem 1.1. The coupled spin-oscillator is a semitoric integrable system, with one single focus-focus 
singularity at m = (0, 0, 1, 0, 0) € S'^ x R^. The semiglobal dynamics around m may be described 
as follows: the linear deviation from exhibiting logarithmic behavior in a saturated neighborhood of m is 
given by the linear map L : M^ — )• M with expression L{X, Y) = ^ X + 5 In 2 Y. In other words, we have 
an equality {S{X, Y))°^ = L{X, Y) + 0{X,Y)^, where {S{X, Y))°° denotes the Taylor series invariant 
at the focus-focus singularity. 

As far as we know, this theorem gives the first rigorous estimate in the literature of the logarithmic 
deviation, and hence the first explicit quantization of the symplectic dynamics around the singularity; we 
prove it in Section 2. The proof is computational but rather subtle, and it combines a number of theorems 
from integrable systems and semiclassical analysis. The method of proof of Theorem II . II ( given in several 
steps) provides a fairly general algorithm to implement in the case of other semitoric integrable systems. 
Moreover, it seems plausible to expect that the techniques we introduce generalize to compute higher order 
approximations, but not immediately - indeed, the linear approximation relies on various semiclassical 
formulas that are not readily available for higher order approximations. In this paper we will also find the 
other invariants that characterize the coupled spin-oscillator (Section 3): the polygon and height invariants; 
these are easier to find. 



Spectral theory for quantum coupled spin-oscillators 

Sections 4, 5 of this paper are devoted to the spectral theory of quantum coupled spin-oscillators. The 
following theorem describes the quantum spin-oscillator. For any h > such that 2 = h{n + 1), for some 
non-negative integer n G N, let Ti denote the standai^d n+ 1-dimensional Hilbert space quantizing the sphere 
S"^ (see Section HUl. 



Figure 1.2: Semiclassical joint spectrum of J, H. We will explain this figure in more detail in Section 4. 



Theorem 1.2. Let S'^ x M? be the coupled spin-oscillator, and (as above) let J, H: M 
commuting smooth functions that define it. The unbounded operators J := Id ® 
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and H = \{x ® u + y ® {^-§^^) on the Hilbert space Ti ® L^(M) C 1?{M?) ® 1? 
commute. The spectrum of J is discrete and consists of eigenvalues in h{^^ + N). 

For a fixed eigenvalue A of J, let E\ := ker( J — Aid) be the eigenspace of the operator J over A. There 
exists a basis B\ of £\ in which H restricted to £\ is given by 
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where < k < n, Iq := ^ + ^, /x := 
The dimension of£\ is /i + 1. 

Finding out how information from quantum completely integrable systems leads to information about 
classical systems is a fascinating "inverse" problem with very few precise results at this time. Section 5 
explains how information of the coupled spin-oscillator, including its linear singularity theory (computed in 
Section 2), may be recovered from the quantum semiclassical spectrum. 

The way in which we recover this linear singularity theory relies on a conjecture for Toeplitz operators, 
which has been proven for pseudodifferential operators. We explain in detail how to do this and formulate 
the following conjecture about semitoric integrable systems: that a semitoric system is determined up to 
symplectic equivalence by its semiclassical joint spectrum, i.e. the set of points in M^ where on the x-axis 
we have the eigenvalues of J, and on the vertical axis the eigenvalues of H restricted to the A-eigenspace 
of J. From any such spectrum one can construct explicitly the associated semitoric system. We give strong 
evidence of this conjecture for the coupled spin oscillators. 
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2 Singularity theory for coupled spin-oscillators 

This section considers semiglobal properties. It is independent of Section 3 which concerns global proper- 
ties. The main goal of this section is to prove Theorem ll.il 

Let (A/, uj, F := (J, H))bea. semitoric integrable system. Recall that a singular point, or a singularity, 
is a point p G M such that rank(dF) (p) < 2, where F := (J, H): Af — )• M^. A singular fiber of the system 
is a fiber of F : A/ — >■ M^ that contains some singular point. 

Let m be a focus-focus singular point m. Let B := F{M). Let c = F{m). The set of regular values of 
F is lnt{B) \ {c}, the boundary of B consists of all images of elliptic singularities, and the fibers of F are 
connected (see ll25l ). 

We assume that the critical fiber Tm := F~^{c) contains only one critical point m, which according to 
Zung [i28il is a generic condition, and let T denote the associated singular- foliation. 

By Eliasson's theorem ifTOl there exist symplectic coordinates (xi, 2:2, ^i, ^2) in a neighborhood U 
around m in which (gi, 52 )> given by 

qi = a;i6 - 2:26, 92 = a;i6 + 2:26, (2-1) 

is a momentum map for the foliation T (in the sense that for some local diffeomorphism q = g o F, so 
the maps q and F have the same fibers); here the critical point m con^esponds to coordinates (0, 0, 0, 0). 
Because of the uniqueness of the S^-action one may chose Eliasson's coordinates ll22l such that qi = J. 

2.1 Construction of the singularity invariant at a focus-focus singularity 

Fix A' G Tm n ([/ \ {m}) and let E denote a small 2-dimensional surface transversal to T at the point A', 
and let Q be the open neighborhood of Tm which consists of the leaves which intersect the surface S. 

Since the Liouville foliation in a small neighborhood of S is regular for both F and q = {qi, q2), there is 
a local diffeomorphism ip of M^ such that q = ipo F, and we can define a global momentum map ^ = ipo F 
for the foliation, which agrees with q on U. Write <1> := {Hi, H2) and Ac := <l>~^(c). For simplicity we 
write (^ = q. Note that Aq = Fm- It follows from (12.11 ) that near m the //i -orbits must be periodic of 
primitive period 27r. 

Suppose that ^ G Ac for some regular value c. Let T2{c) > be the time it takes the Hamiltonian flow 
associated with H2 leaving from A to meet the Hamiltonian flow associated with Hi which passes through 
A, and let ri(c) G IR/27rZ the time that it takes to go from this intersection point back to A, hence closing 
the trajectory. We denote by 7c the corresponding loop in Ac. 

Write c = (ci, C2) = ci + ic2, and let In z for a fixed determination of the logarithmic function on the 
complex plane. Let 

0-1 (c) = ri(c)-9(lnc) 
0-2 (c) = r2(c) + K(lnc), 




S\A) 
Figure 2.1: Singular foliation near the leaf J>„, where -S'^(^) denotes the 5'^-orbit generated by }i\ = J. 



where K and 9 respectively stand for the real an imaginary parts of a complex number. Vu Ngoc proved 
in |[23l Prop. 3.1] that a\ and 02 extend to smooth and single- valued functions in a neighbourhood of and 
that the differential 1-form 

a := CTi dci + (T2 dc2 

is closed. Notice that if follows from the smoothness of a2 that one may choose the lift of T2 to M such that 
(72(0) G [0, 2tt). This is the convention used throughout. Following ll23l Def. 3.1] , let S be the unique 
smooth function defined around G M^ such that 

dS = a, 5(0) = 0. 

The Taylor expansion of S at (0, 0) is denoted by (5)°°. 

The Taylor expansion (5)°° is a formal power series in two variables with vanishing constant term, and 
we say that {S)°° is the Taylor series invariant of {M, ui, (J, H)) at the focus-focus point c. 

2.2 The coupled spin-oscillators 

Let S"^ be the unit sphere in M?' with coordinates (x, y, z), and let M^ be equipped with coordinates (n, v). 
Recall from the introduction that the coupled-spin oscillator model is the product S"^ x M^ equipped with the 
product symplectic structure ugi © cjq given by dO f\(iz®duf\dv, and with the smooth Poisson commuting 
maps J, //: M — )• M given by J := [u^ + v^)/2 + z and H := ^ {ux + vy). Sometimes we denote 
the coupled spin-oscillator by the triple {S'^ x M^, ujg2 ® wq, (./, H)). A simple verification leads to the 
following observation. 

Proposition 2.1. The coupled spin-oscillator (S^ x R^, 00^2 © ojqAJ, H)) is a completely integrable 
system, meaning that the Poisson bracket {J, H} vanishes everywher^ 

In addition, the map J is the momentum map for the Hamiltonian circle action of S^ on S'^ x M? that 
rotates simultaneously horizontally about the vertical axes on S'^, and about the origin on M?. 

The singularities of the coupled spin— oscillator are non-degenerate and of elliptic-elliptic, transversally- 
elliptic or focus-focus type. It has exactly one focus-focus singularity at the "North Pole" {{0, 0, 1), (0, 0)) G 
5^ X M? and one elliptic-elliptic singularity at the "South Pole" ((0, 0, —1), (0, 0)). 



equivalently the Hamiltonian vector field Xj is constant along the flow of Xh 



Corollary 2.2. The coupled spin-oscillator (5^ x M^, ujg2 © wq, {J, H)) is a semitoric integrable system. 

Computing the Taylor series invariant at the focus-focus singularity is rather involved. At this point we 
are able to compute the first two terms ai, 02 (for the coupled spin-oscillators). Even in this case one has to 
do a delicate coordinate analysis of flows involving Eliasson's coordinates, and the computation of various 
integrals. 
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Figure 2.2: Spin model with momentum map z. Here (6*, z) are the angle-height coordinates on the unit 
sphere 5^. 



2.3 Set up for coupled spin-oscillators — Integral formulas for singularity invariant 

Throughout we let M = S*^ x M^ and F = (J, H). In this set up stage we introduce the 1-forms ki c and 
k;2,c in terms of which the Taylor series in defined in Il23l . and we recall limit integral formulas for the Taylor 
series invariant. Then we introduce the limit theorem proved in the semiclassical paper |[22l Proposition 6.8], 
which will be the key ingredient for the computation . 

The formulas that we present here do not correspond to the exact statements in the corresponding papers, 
but can be immediately deduced from it assuming the context of the present paper. 



The one forms ki c and K2,c- As usual, we denote by Xq^ the Hamiltonian vector field generated by qi, 
i = 1, 2. Let c be a fixed regular value of F. Let ki c G ri^(Ac), ki^c G f^^(Ac) be the smooth 1-forms on 
the fiber Ac := F~^{c) corresponding to the value c defined by the conditions 



l^l,c{Xqi) '■— —1, l^l,c{'^c 



92 y 



0, 



and 



K2,c('Vgi) := 0, K2,ciXq^) ■= -1- 



(2.2) 



(2.3) 



Note that the conditions in (12.21 ) and (12.31 ) are enough to determine ki_c and K2,c on Ac because Xg-^ , Xg^ 
form a basis of each tangent space. 

We will call ni^, K2,o the corresponding form defined in the same way as ki c, K2,c, but only on Ao\{?n}, 
where m = (0, 0, 1, 0, 0) is the singular point of the focus-focus singular fiber Aq. 

Remark 2.3 The forms ki c, K2 c > * = 1, 2 are closed. See also |[22l Section 3.2.1]. 



Limit integral formula for Taylor invariants. The following result will be key for our purposes in the 
present paper. 

Lemma 2.4. Let (S) G IK[[X, Y]] be the Taylor series invariant of the coupled-spin oscillator. Then the 
first terms of the Taylor series are given by the limits of integrals ai = limc-s>o ( / i^i.c + arg(c) I and 

02 = limc^o f Xy^ K2,c +ln |cM. 

Proof. It follows from the definition of the dynamical invariants ti(c) and T2(c) in Section |2] and the defi- 
nition of Ki c and K2,c in (12.21 ) and (12.31) respectively that Ti{c) = f Ki^c, i = 1, 2. The first two terms of 
the Taylor series invariant ai (0) and 0-2(0) where ui = ti + arg(c) and cr2 = T2 — In |c|. 

Since o"i and (T2 are smooth, we have that oi = (Ti(0) = limc_j.o ( / ^i,c + arg(c) ) and 02 = (72(0) = 

Localization on the critical fiber. On the other hand, we have the following [22^, Proposition 6.8] result 
proved by the second author. 



Theorem 2.5 ( 11221 ). Let 70 be a radial simple loop. The integrals in Lemma \2A\ are respectively equal to 

/^i.c + arg(c) = lim / ki + (t^ - 6'b) , (2.4) 



and 



|.Bo(t):=7o(l-t) 

02 := lim ( / K2,c+ln|c|)= lim ( / ^2,0 +ln(r^oPBo) ) ' (2-5) 

' h. ' (s,t)^(0,0) \ } A<i:=-y^{s) 



C-5>0 



where for any point A in M close to m with Eliasson coordinates (xi, X2, ^i, ^2) ^s defined in equation 
d2.iD . we denote by (r^, tA, PA, Oa) the polar symplectic coordinateso of A, i.e. {va, tA) are polar coor- 



dinates corresponding to (xi, X2) and {pA, Oa) are polar coordinates corresponding to (^1, ^2)- 

2.4 Computation of integral limit formulas for coupled spin-oscillators 

Now, in order to apply Theorem l2.5l we need to find the curve 70, as well as the 1-form k and the coordinates 
(r, 6, p, a), both of which are defined on Aq. First we describe a parametrization of Aq, and then we use 
this pai^ametrization to define 79. We have divided the computation into five steps. 

Stage 1 — Eliasson's coordinates (xi, X2, fi, ^2) 

We find explicitly symplectic coordinates (xi, £2, £,1, C2) ^ M = S'^ x M^ in which the "momentum map" 
F: M — )• M^ for the coupled spin-oscillator has the form (12. lb . up to a third order approximation, i.e. up to 
(0(xi, X2, ii, 6))^. For brevity write 0(3) = (0(xi, xa, 6, 6))^- 



^These coordinates (ta, tA, Pa, Oa) should not be confused with the coordinates (r, t, p, 9) without the subscript, which are 
coordinates in R^ x S^. 



Lemma 2.6. Consider the map (j) '■ T(q^ q^ q, o) ^^ ~^ '~^{o, o, i, o, o) i^'^ ^ ^^^) given by 

4>{xi, X2, li, I2) = {v ■■= ^(^2 + li), a: := ^(^2 - li), w := 7^(~^i + ^2), V ■= ^(^1 + ^2))- 

TIw map (j) is a linear symplectomorphism, i.e. an automorphism such that (ffVt = ujq, where wq = dxi A 
d^i © dx2 A d^2 is the standard symplectic form on R^, and i7 = (a'52 © du A dv) |^y (s'^xM?) (recall 

01^2 is the standard symplectic form on S'^). In addition, (p satisfies the equation Hess(-F) o (j) = [qi, (72). 
where F := B o {F - F{m)) = B o (F - [1, 0)) : M -)■ M?Jor the matrix ^ := ( J \ 

In the above statement, we identify a Hessian with its associated quadratic form on the tangent space. 

Stage 2 — Curve and Singular Fiber Parametrization 

Parametrization of Aq. Let's now parametrize the singular fiber Aq := F~^(l, 0), where F = (J, H) as 
usual. This singular fiber Aq corresponds to the system of equations J = 1 and H = 0, which explicitly is 
given by system of two nonlinear equations J = [u"^ + v^)/2 + z = and H = ^{ux + vy) = 0. on the 
coordinates {x, y, z, u, v) on the coupled spin oscillator M = S"^ x R^. 

In order to solve this system of equations we introduce polar coordinates u + iv = re'* and x+iy = pe'^ 
where recall that the 2-sphere S'^ C M.^ is equipped with coordinates (.x, y, z), and M^ is equipped with 
coordinates (n, v). 

For e = ±1, we consider the mapping S^ : [— 1, 1] x M/27rZ — ;■ M^ x 5^ given by the formula 
S,{p) = {r{p) e<P\ {p{p) ei^(p)'^(p))) where p = (z, 9) G [-1,1] x [0,2^) and 

'r(p) = V2(l-5) 
t{p) = ~e + e^ 

p{p) = Vi -52 
e{p) = 9 

^z{p) = z. 

Proposition 2.7. The map S^, where e = ±1, is continuous and S^ restricted to (—1, 1) x M/27rZ is a 
diffeomorphism onto its image. If we let Aq := 5e([— 1, 1] x IR/27rZ), then Aj U Aq = Aq and 

Ai n A2 = ({(0, 0)} X {(1, 0, 0)}) U {C2 X {(0, 0, -1)}) , 

where C2 denotes the circle of radius 2 centered at (0, 0) in M^. Moreover, S^ restricted to (—1, 1) x ]R/27rZ 
is a smooth Lagrangian embedding into M^ x S"^. 

Proof. On the one hand we have that z"^ = 1 — x"^ — y"^ = I — p^. The expressions for the maps J and H 
in the new coordinates (r, t, p, 9) are 

J=\r^^ Vl^, H=^-^cosit-9). (2.6) 



In virtue of the formula for H in the right hand-side of (12.61 ). if H = then r = or p = or t — = 
(modyr), which leads to thi^ee separate cases. The first case is when r = 0; then J = ±^1 — p"^ = 1, and 



hence p = 0. Hence the only solution is (u, v, x, y, z) = (0, 0, 0, 0, 1). The second case is when p = 0; 
then either z = 1 and r = 0, or z = — 1 and r = 2. Hence the set of solutions consists of (0, 0, 0, 0, 1) and 
the circle r = 2, p = and z = —1. Finally, the third case is when t — 6 = -Kmod vr); because J = 1 and 
H = 0,it follows from the formula for z above and the left hand-side of ( I2.6I ) that r^ = 2{1 — z). Hence the 
set of solutions Aq is equal to the set of points (re'*, p e'^) such that 

= ^/2{l-z), z G [-1, 1] 

= t-f or e = t + ^, tG[0, 27r) (2.7) 

lp= \/l - z^ 



This case contains the previous two cases, which proves statement (3) part (i) in virtue of expression (I2.5I ). 
The other statements are left to the reader. D 

Remark 2.8 The singular fiber Aq consists of two sheets glued along a point and a circle; topologically Aq 
is a pinched torus, i.e. a 2-dimensional torus S^ x S^ in which one circle {p} x S^ is contracted to a point 
(which is of course not a a smooth manifold at the point which comes from the contracting circle). 



The radial vector field Xy, on Aq. 

Proposition 2.9. Let Xg. be the Hamiltonian vector field of qi (which recall is defined in saturated neighbor- 
hood of the singular fiber AqJ. On the singular fiber Aq, the vector fields Xq-^ , Xj and Xq^ , Xh are linearly 
independent, precisely: Xq^ = Xj, Xq^ = 2 Xh- In particular the vector field Xh is radial. 

Proof. It follows from Eliasson's theorem that there exists a smooth function h such that q = h o F and 
d/i(0) is the invertible 2 by 2 matrix B in Lemma l2!6l 
Then on Aq we have that 

y=^Xj + ^XH, i = l,2. (2.8) 

Because the coefficients ai^e constant along Aq, it is sufficient to do the computation at the origin. At 
the origin the computation is given by the matrix B in Lemma |2?6l so we have that ^j{0) =, ^^^(O) = 
0, ^(0) = and |^(0) = 2. The proposition follows from dH]). D 

In the following section we will need to use explicitly the Hamiltonian vector field Xh, and therein it 
will be most useful to a have the following explicit coordinate expression. 

Lemma 2.10. The Hamiltonian vector field Xh of H is of the form 

yd X d —yu -\- xv d z{xu + yv) d 



2du 2dv 2 dz 2(1 - z2) 06' 

Proof. For this computation let us use coordinates [u, v, z, 6) as a parametrization of M^ x S'^. 

The coordinate expression for the Hamiltonian H is H = ^{xu + yv) = ^{p cos 6u + psinOv), Then 
the Hamiltonian vector field Xh is of the form Xh = a-^ + b-^ + c-^ + d-^, where since the symplectic 
form on M^ x S^ in these coordinates isduAdt;+d0Adz, the function coefficient a (which will be important 
later in the proof) is given by 
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and the other function coefficients are given by 6 = — ^ = pcos{9) = — |, c = ^ = f^" sm{6)u + 
cos{9)v) = ^^^ and d = -||. 

We need to compute d expliciriy. Since || = because the angle 9 does not depend on the height z, 
and 4^ = — 7^=^, we have that 

dx dx dp dx 89 dx dp —xz 

— = L- -I = L. = (2.10) 

dz dp dz d9 dz dp dz p^ 



+ ^^ = ^^ = ^ (2-11) 



dy dy dp dy d9 dy dp —yz 
dz dp dz d9 dz dp dz p 

It follows that from ( I2.10I ) and ( I2.11I ) that the function coefficient d is given by 



dH dH dx dH dy u —xz v —yz z[xu + yv) z{xu + yv) 

~ ~~dz ~ ~~d^'dl^~dy"d^ ~ 2 ^ "^ 2 V" ~ V ~ 2(1-^2) 



D 



Definition of a simple "radial" loop in Aq. In order to apply the theorem it is enough to take 70 to be an 
integral curve of the radial vector field Xh- 

We define 70 as the simple loop obtained through the parametrizations S+ and S^ by letting z run from 
—1 to 1 and back to —1, respectively. For instance, one can use the formula 

. ^ / 5i(-l + 4s, -f) ifO<s<i; 



52(3-4s, f) ifi<s<l. 
Corollary 2.11. Along the curve 70 we have 

yd yu d 



Xu 



r, ■ (2.12) 

70 2 du 2 dz 



Proof. We use the notation of Lemma [2.101 Along 70 we have f = 0, x = and 9 = tt or 9 = ^. Hence 
a = |,6 = 0, ,c = — ^,d = 0. Therefore the vector field Xh along the curve 70 is given by (12.121 ). D 

Using Corollary 12.111 we describe the very explicit relation between the curve 70 and the Hamiltonian 
vector field Xh- 

Proposition 2.12. The curve 70 : [0, 1] — s- M is an integral curve of Xh- 

Proof. Since by construction the vector field 5'*(^) is tangent to the curve 70, it is enough to show that 
S*{-^) is colinear to Xh ai^e colinear at each point. 
A computation gives that 

s (^\ = .^- ^ I- + ^ I. a 13) 

*\dz) dz ^2(1 -z) dr^ ^Y^£ldp ^- ^ 

On the other hand 

« = V2(l-z), (2.14) 
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and since (r, t) aie polar coordinates for {u, v), -^ = cos t ^ +sin t ^, which at t = gives that ^ = ^■ 
Therefore, because at t = the last factor of ( I2.13I ) is zero, we conclude from ( I2.14I ) that 

\dzJ dz u du 

It follows from (12.12b that Xh = ~^ '^*[m)' which shows that Xh and S* ( ^ ) are cohnear at every 
point, as desired. 

D 

Stage 3 — Integration in linearized Eliasson's coordinates 

Let be a local symplectic map such that goFocl) = qon M^, as given by Eliasson's normal form theorem. 
The integrals in (..) are defined in terms of the corresponding canonical coordinates (xi, X2, Ci, ^2) in I^^- 

Because our computation is local, we can use instead the linearized coordinates that we have de- 
fined in Lemma 12.61 More precisely, one can always choose </> such that the tangent map d(o, 0, 0, 0) </* • 
T(o,o,o,o) J^^ — ;• T(o,o,i,o,o) '^'^ ^ ^"^ is equal to 0, and this gives local coordinates (xi,X2,^i,^2) in a 
neighborhood of m, such that B o F{xi, £2,^1,^2) = '?(a^i,^2,'?i,'?2) + 0{3). 

Note that these coordinates are not symplectic, except at m. 



Lemma 2.13. The integral f l2.5D gives us the same result when computed in linearized coordinates, i.e. upon 
replacing va by ta, ^a by Ia, Pa by pA and 9a by 6a- 

Proof. Since r\ = x\ + x\, then 

r\ = xl + xl = xl + xl + 0(3) =r\ + 0(3) (2.16) 

We know that -§^ = 0(1), and therefore it follows from (IZT6l ) that 

ln(fi)=ln(r^ + 0(3))=ln(l + ^)+In(r^)=ln(l + 0(l))+ln(ri) = 0(l)+ln(ri).(2.17) 

Similarly ln(p|) = 0(1) +ln(p|). Hence In(rAPB) =ln(rA) +ln(pB) =\n{fA) +Mpb) =\i'i{fAPB) + 
0(l).Then 

. ^™}. „ In(^AoPBo) -ln(rAo/5B„) = 0. (2.18) 

(s,t)->(0, 0) 

It follows from expressions (12.51) and (12.181 ) that 

/>Bo=7o(1-sb) s 

02= lim ( / «:2,o+ln|fAoPBol)- (2-19) 

(sa,sb)^(0,0) \JA(^=^o{sa) ' 

This concludes the proof. D 
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Stage 4 — Computation of the first order Taylor series invariants ai and 02 

In order to compute the integrals in ( I2.19I ) we can replace 70 by any integral curve of X^ with the same 
endpoints. Thus, let 7 be a solution to 7 = Xh o 7. By definition, for any 1-form k, 

-B():=7(s2) rS2 rS2 

^= I S{s)(7(s))ds= / K^i^s){^H{l{s)))ds. (2.20) 

^o:=7(si), along7 J si J si 

Theorem 2.14. Let [S] G IR[[^, ^]] be the Taylor series invariant of the couple-spin oscillator Then the 
first coefficient of the first term of the series is given by ai = ^. The second coefficient of the first term of 
the first order Taylor series invariant is 02 = 5 In 2. 

Proof We have divided the computation of 02 in several steps. 

Step 1: Set-up of the integral 0/^2,0- We need to compute expression (12.191 ). 



Let a be given by (12.91 ). 

In view of (12.121 ). the path 7 between Aq and Bq can be parametrized by the variable u. This means 
that the path 7 is obtained by first increasing it up to u = 2 on the first sheet (parametrized by ^i) and then 
decreasing u on the second sheet (parametrized by S2). 



By LemmaEHwe know that Xg^ = 2Xh and hence {K2,o)-y{s){'^Hi'y{s))) = ^'°^^W^^ ''a^'^^ '" _ By 

-Bo 

Aq , along 7 ' 



definition of ^2,0 we know that ^2,0 (^92) = ~1 ^^d hence it follows from (12.201 ) that /^"^lon ^2,0 



J"*^ ^. Since ^ is equal to a = | we have that 



^0 nds f^ du r^ du 

AJ2,o = / ^ = / —f^ + / ^r^, (2.21) 

Ao,along7 Jsi 2 J^^y+{u) J2 y~[U) 

where y±{u) is the y-coordinate along the part of the curve 70 which corresponds to the parametrization 
S±, respectively. Our next goal is to compute expression (12.211) . 

Step 2: Computation of expression ii2.21\) . Now, y = psm{6) = ±p. 

Now let us express the dependence of y in u along the path 7. By the equation J = i(u^+t;^) + z = — 1, 

2 

which is always true along the singular fiber, we have that, since i; = 0, ^ + z = 1, or in other words, 

2 

z = 1 — ^.It follows from this equation that 



n^ 



y^ = ±p = ±\/l-z2 = ±Wi _ (1 _ )2 = ±uJl - — since -u > 0. (2.22) 

On the other hand, note that the function G(t) = In ( -^^ + tan i ) is a primitive of the function g{f) = 

■^^. Then by equation (12. 221) .using the change of variable u/2 = cos t, and then applying the fundamental 
theorem of calculus we obtain^ 

'■^ du f^ du - , ^ 

m I h tan t 

V cos t 



ui y+ Jui ^^/i ^ 



2 








^2 2 


1 


?f2xi 


2 


^ — 


In 


- + -\ 


1 1 - 


-^) 




h 




\u u \ 




4 J\ 


Ml 



'The integral is equal to when u — 2 
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and simplifying this expression we then obtain 



:'^^,nl^)+lJl + f^). (223) 

„^ y+ \uiJ V V 4 7 



The goal of this proof is to compute ai , which by (12.51 ) is equal to the limit 

lim / K2fi +lii{rAoPBo] 

and precisely because this limit exists, we may calculate it along the diagonal values given by u = ni = U2. 
Then it follows from equation ( I2.23I ) that 



Ao Ju,y+ J2 y~ Ju y+ V Vn/ V V aJJ 

This concludes this step. 

Step 3: Computation of the logarithm factor InlrAoPBo) ■ 

From the notation of Stage 1 we have that f\ = x\ + £2 ^^i^ that p\ = ii + ^2- Using Lemma [Z6l we 
find that f\ = \{x^ + y^ + v? + v'^) + {—uy + vx) and p\ = \{x'^ + y"^ + u^ + v"^) + {uy — vx). 

We need to compute taq and pbq ■ The points Aq and Bq are in the path 70 and Aq : = {uao , vaq , Oao j ^Ao ) 

V? V? 

{uao, 0, f , 1 - -^), and Bo := {ub^, vbq, Obq, zbo) = (^Ao, 0, ^, 1- -^). 

With this information we can compute j-Aq and p^y using expression (12.221) and recalling that x = u = 
along 7: 



fl, = \(-'-"^ + -')-'f^ = f(2-^-2/lT^). (2.25) 

where here we have also used p^ = l — z^ = l — (1 — ^)^ = n^ — ^. And we also have that 

/5|o=^io- (2-26) 

It follows from (IIl5l) and (II!26l ) that 

HtaJb,) = \Hr\A) = \^<^m) =l^(^io) =1" (y(2 - X + 2^1 -y)) 



and therefore that 



2 

u , , , u 



ln(fAoM,) = 21n( — ) +ln(2 - _ + 2^1 - -)• (2.27) 



This concludes the computation of the logarithmic factor. 
Step 4: Conclusion. It follows from (|23] ). (12.241) and (12.271 ) that 

02 = lim ( / K2,o +ln(fAo/5B(,) 



lim((21n(-) + 21n(l + ^l--) + 21n(^)+ln(2-- + 2^1-- 

21n2 + 21n2-ln2 + 21n2 = 51n2. (2.28) 
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So we have proven that 02 = 51n2 as we wanted to show. 

In order to find ai, note that the following hold: n > 0, v = 0, 6 = ^ or ^, p = \/l — z"^, z 

1 - ^, p = Jv? - ^. In this case xi = ^, X2 = ^, and therefore ^ = f . Similarly ^1 

-\^-, C2 = -^^ = — '^1) and hence a = ^. It follows that Oaq — a_Bo — f • Therefore by Theorem 

Bo=7o(i-t) , ^ ^ 



ai = lim / Ki + {Oa - olb) . 

(s,t)^{0,0) V 7^0=70(1) ^ 2 

Here we are using that because kq{Xh) = and 70 is tangent everywhere to Xh so one has that 

lim ( / kq) = 0. 

is,t)-,iO,0)\JAo='yo{l) ' 

(See also the paragraphs before Theorem 12. 14| ). This concludes the proof. D 

Theorem 11.11 follows from Theorem 12. 141 

Remark 2.15 It is plausible that our proof technique generalizes to compute the higher order terms of 
the Taylor series invariant, but not immediately, as we rely on the limit theorem proved in ll22]| which only 
applies to the first two terms. The computation provides more evidence of the fact that from a dynamical 
and geometric view-point focus-focus singularities contain a large amount of information. 



3 Convexity theory for coupled spin-oscillators 

The plane M^ is equipped with its standard affine structure with origin at (0,0), and orientation. Let 
Aff(2,]R2) :=GL(2,R2)xM2be the group of affine transformations ofM^. LetAff(2,Z) := GL(2,Z) xM^ 
be the subgroup of integral-qffine transformations. 

Let T be the subgroup of AflF(2, Z) of those transformations which leave a vertical line invariant, or 
equivalently, an element of T is a vertical translation composed with a matrix T^, where A; G Z and 



^'-{i ;)"°^<'' 



Let io C M^ be a vertical line in the plane, not necessarily through the origin, which splits it into two half- 
-spaces, and let n G Z. Fix an origin in i. Let t" : M^ — ;• M^ be the identity on the left half-space, and T" 
on the right half-space. By definition t" is piecewise affine. A convex polygonal set A is the intersection 
in M? of (finitely or infinitely many) closed half-planes such that on each compact subset of the intersection 
there is at most a finite number of corner points. We say that A is rational if each edge is directed along a 
vector with rational coefficients. For brevity, in this paper we usually write "polygon " instead of "convex 
polygonal set" . 

3.1 Construction of the semitoric polygon invariant 

Let £ be a vertical line through the focus-focus value c. Let B^ := 1tA{B) \ {c}, which is precisely the set of 
regular values of F. Given a sign e G { — 1, +1}, let ^"^ C ^ be the vertical half line starting at c at extending 
in the direction of e : upwards if e = 1, downwards if e = — 1. 
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In Th. 3.8 in Il25l it was shown that fore G {— 1,+1} there exists a homeomoiphism f = f^: B ^M?, 
modulo a left composition by a transformation in T, such that f\(^B\e') is a diffeomorphism into its image 
A := f{B), which is a rational convex polygon, flmAi^) i^ affine (it sends the integral affine structure of 
By. to the standai^d structure of M^) and / preserves J: i.e. f{x, y) = {x, f'^'^\x, y)). f satisfies further 
properties ifTSl . which are relevant for the uniqueness theorem proof. In order to arrive at A one cuts 
(J, H){M) C M^ along the vertical half-lines i,'^. Then the resulting image becomes simply connected and 
thus there exists a global 2-torus action on the preimage of this set. The polygon A is just the closure of the 
image of a toric momentum map corresponding to this torus action. 



S-^x 



{.IH) 




(4,ei = +l) 



Figure 3.1: The coupled spin-oscillator example. The middle figure shows the image of the initial moment 



map F = (J, H). Its boundary is the parametrized curve (j(s) 



.h{s) 



s G [1,00). 



2s '"-V-^/ — -^^7572^ 

The image is the connected component of the origin. The system is a simple semitoric system with one 
focus-focus point whose image is (1, 0). The invariants ai^e depicted on the right hand-side. The class of 
generalized polygons for this system consists of two polygons. 



We can see that this polygon is not unique. The choice of the "cut direction" is encoded in the signs 
e, and there remains some freedom for choosing the toric momentum map. Precisely, the choices and the 
corresponding homeomorphisms / are the following : 



(a) an initial set of action variables /o of the form (J, K) near a regular Liouville torus in 1251 Step 2, 
pf. of Th. 3.8]. If we choose /i instead of /o, we get a polygon A' obtained by left composition with 
an element of T. Similarly, if we choose /i instead of /o, we obtain / composed on the left with an 
element of T; 

(b) an integer e G {1, —1}. If we choose e' instead of e we get A' = t„(A) with u 
Prop. 4.1, expr. (11)]. Similarly instead of / we obtain f = t^o f. 



(e - e')/2, by ||25 



Once /o and e have been fixed as in (a) and (b), respectively, then there exists a unique toric momentum 
map fj, on Mr := F~^{lntB \ ie) which preserves the foliation T, and coincides with foo F where they are 
both defined. Then, necessarily, the first component of fi is J, and we have f-i{Mr) = A. 
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We need now for our purposes to formalize choices (a) and (b) in a single geometric object. Let 
Polyg(M^) be the space of rational convex polygons in M^. Let Vert(]R^) be the set of vertical lines in 
M^. A weighted polygon (of complexity 1) is a triple of the form A^ = [A, ix, ej where A € Polyg(IR^), 

£ G Vert(IR^), and e E { — 1, 1}- Let G := {—1, +1}- Obviously, the group T sends a rational convex 
polygon to a rational convex polygon. It corresponds to the transformation described in (a). On the other 
hand, the transformation described in (b) can be encoded by the group G acting on the triple A^ by the 
formula 

e'-(A,£x,e) = (tu{A),ix,e'e] 



where n = (e — e')/2. This, however, does not always preserve the convexity of A, as is easily seen when 
A is the unit square centered at the origin and Ai = 0. However, when A comes from the construction 
described above for a semitoric system (J, H), the convexity is preserved. Thus, we say that a weighted 
polygon is admissible when the G-action preserves convexity. We denote by >VPolyg(M^) the space of all 
admissible weighted polygons (of complexity 1). The set G x T is an abelian group, with the natural product 
action. The action of G x Ton yVPolyg(M^), is given by: 



(6', t).(a, £a, e) = (i«MA)), 4, e'e), 



where u = {e — e')/2. We call a semitoric polygon the equivalence class of an admissible weighted polygon 
under the {G x T) -action. 

Let A be a rational convex polygon obtained from the momentum image (J, H){M) according to the 
above construction of cutting along the vertical half-line £'^. 

Definition 3.1 The semitoric polygon invariant of (M, uj, (J, H)) is the semitoric polygon equal to the 
(G X r)-orbit (G X T) • f A, I, e) G WPolyg(M2)/(G x T). 



3.2 The semitoric polygon invariant of coupled spin-oscillators 

Proposition 3.2. The semitoric polygon invariant of the coupled spin-oscillator is the (G x T)-orbit con- 
sisting of the two convex polygons depicted on the right hand-side ofFieure \3.1\ 



Proof. As shown in Figure 13.11 a representative of the semitoric polygon invariant is a polygon in M? with 
exactly two vertices at (—1, 0) and (1, 0), and from these two points leave straight lines with slope 1 (the 
other possible polygon representative has vertices at (—1, 0) and (1, 2)). One finds this polygon simply by 
combining the information about the isotropy weights at the left corner of the polygon (an elliptic-elliptic 
critical value) [25, Prop. 6.1], together with the formula given in |[25l Thm. 5.3], in which the relation 
between isotropy weights and the slopes of the edges of the polygon is described using the Duistermaat- 
Heckman function. D 

3.3 Classification theory for coupled spin-oscillators 

The authors have recently given a general classification of general semitoric integrable in dimesion 4 lITSl . 
|[T9l in terms of five symplectic invariants; the reader familiar with these works can easily that two of these 
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invariants do not appeal" in the case of coupled spin-oscillators, and we state the uniqueness theorem therein 
in this particular casqj 

Consider a focus-focus critical point m whose image by (J, H) is c, and let A be a rational con- 
vex polygon con^esponding to the system (M, w, (J, H)). If /i is a toric momentum map for the system 
(Af, uj, (J, H)) corresponding to A, then the image ii{'m) is a point in the interior of A, along the line I. 
We proved in fTS^ that the vertical distance h := n{m) — min^g^nA t^2{s) > is independent of the choice 
of momentum map fi. Here 7r2 : M^ — >■ M is 7r2(ci, C2) = C2- 

Theorem 3.3 (consequence of Th. 6.2, HU). Let (Af, uj, (J, H)) be a A-dimensional semitoric integrable 
system with exactly one focus-focus singularity. The list of invariants of (Af, w, (J, H)) consists of the 
following items: (i) the Taylor series invariant {S)°° at the focus-focus singularity m; (ii) the semitoric 
polygon invariant; (Hi) the volume invariant, i.e. the height h > ofm. Two 4-dimensional simple semitoric 
integrable systems (Afi, wi, (Ji, Hi)) and {M2, c<J2, (<f2) -^^2)) with exactly one focus-focus singularity are 
isomorphic if and only if the list of invariants (i)-(iii) of ( Af 1 , wi , ( Ji , Hi ) ) is equal to the list of invariants 

(i)-(iii)of{M2,0J2,{J2,H2)). 

Theorem 3.4. The coupled spin-oscillator has the following symplectic invariants: (i) first terms of the 
Taylor series invariant: ai = ^ and a2 = 5 In 2; (ii) semitoric polygon invariant: [G y. T) ■ A^, where 
Aw is either the upper or lower weighted polygon depicted on the right-most side ofFieure \3.1[ (Hi) volume 
invariant: h = 1. 

Proof. The semitoric polygon invariant and the first terms of the Taylor series invariant were computed 
previously. The height of the focus-focus point of the system in the polygon is equal to half of the Liouville 
volume of the submanifold of Af given by the equation J = 1. This is because the functions H and 
J are symmetric about the J-axis of M^ in the sense that J(x, y, z, u, v) = J{x, y, z, —u, —v) and 
H{x, y, z, u, v) = —H{x, y, z, —u, —v). Here there is no need to compute anything because the volume 
of the submanifold given by J = 1 in Af is just the length of the vertical slice of the polygon at J = 1, 
which is 2, and hence the height of the focus-focus point of the system is hi = 1, and the image of the 
focus-focus point in the polygon is (1, 1). D 

4 Spectral theory for quantum spin-oscillators 

In this section, we use the notation of the previous sections J = " ^^ + z and H = ^{xu + vy). Our goal 
in this section is to quantize this example and analyze its semiclassical spectrum. 

First we quickly review the process of assigning a quantum system to a classical system. Loosely 
speaking, a quantum integrable system is a collection of commuting self-adjoint operators on a Hilbert 
space. Quantization is a process that takes a classical phase space (here, a symplectic manifold Af ) to a 
Hilbert space Af, and classical Hamiltonians / G C°°(Af) to self-adjoint operators / acting on Af. The 
quantization of symplectic manifold is often called geometric quantization. See the recent book by Kostant- 
Pelayo |[T5l for a survey. Quantizing Hamiltonians involves more difficulties. For instance, we need the map 



■*The first of these invariants is thie number of focus-focus singularities. The last of these invariants, the so called twisting index 
invariant, is a rather subtle topological invariant which measures how the topology near a focus-focus singular fiber relates to the 
topology near the other focus-focus fibers. Hence the invariant only appears when there is more than one focus-focus singularity, 
and in the following we shall not mention it. The twisting-index expresses the fact that there is, in a neighborhood of any focus- 
focus point Ci, a privileged toric momentum map v. This momentum map, in turn, is due to the existence of a unique hyperbolic 
radial vector field in a neighborhood of the focus-focus fiber. Therefore, one can view the twisting-index as a dynamical invariant. 
This is an important invariant in the general case, see 1181 . 

18 



/ I—)- / to be a Lie algebra homomorphism, at least at first order : if the classical system is given by two 
Poisson commuting functions /, g then the quantum system is given by two operators /, g such that 

f[/, 5]=0 iT^od {0{H)). (4.1) 

Such a quantization is well-knowrlj to exist when M = M^", and more generally on a cotangent bundle 
M =T* X, using ^-pseudodifferential quantization Q. Quantizing compact symplectic manifolds is also 
possible under an integrality condition (the existence of a so-called prequantum line bundle), using Toeplitz 
quantization lH. However, because of the remainder in (14. Il l, it is not known whether a classical integrable 
system can always be quantized to a true quantum integrable system. Very recently, in the algebraic setting, 
the relevant obstruction was defined ifTTI . In the coupled spin-oscillator example, like in many known 
systems, an exact quantization can be found by hand. 

A well-known example is the harmonic oscillator in M^. The harmonic oscillator is given by M = M^ 
with coordinates (u, v) and Hamiltonian function on it N{u, v) = " ^" . The self-adjoint operator N in the 
Hilbert space \? (M) given by iV = — -y -^ + ^ is the standard Weyl quantization of the Hamiltonian N. 
The spectrum of N is discrete and given by {h{n + i) | n G N}. The eigenfunctions are Hermite functions. 
This operator will be used as a quantum building tool in the sequel. 

4.1 Quantization of M} and the Harmonic Oscillator 

We shall view S"^ as a reduced space of M"^ ~ C^ under the coordinate identification zi = xi + i^i, 
Z2 = X2+ 1^2- On M"^ we consider the well-known harmonic oscillator, -^(zi, Z2) = ^ which has a 
27r -periodic flow generating a Hamiltonian 5^-action t ■ {zi, Z2) = (2:1 e~'*, 2:26"'*). 

The space Ye := {L = E}, for any value £^ > 0, is of course the euclidean 3-sphere 5"^/-^ C M^ 

of radius \/2E. It is well known that the reduced space {L = E}/S^ is 2-sphere, and the fibration map 
{L = E} — ;• {L = E}/S^ is the standard Hopf fibration. More precisely, we may represent this reduced 
space as the euclidean sphere 5*1, ,2 '^ ^'^ of radius E/2. Denoting by (x, y, z) the variables in M^, we have 
the following useful formula for the Hopf map, which will be used for quantization : 

X = ^{ziZ2)/2 

y = '^{ziZ2)/2 

z = (|zi|2-|z2|')/4. 



The usual quantization of M^ is the Hilbert space T-L^a = L^(M^). The Weyl quantization of the Hamil- 
tonian function L is the unbounded operator L := — ^ ( j-r + j-j ) + ^^ ^^^ ■ 



The spectrum of L is given by spec(L) = {h{n + 1) | n G N}. To see this, define the operator Lj 



Y I j-7 I + -j; acting on L (R^;^). We can write L = Li + L2. Note that the spectrum of Lj is 

spec(Lj) = {h{nj + -) | Uj G N}. (4.2) 

Therefore we deduce that the spectrum of L is given by {h{ni + 77-2 + 1) | ni G N, n2 G N}, and the formula 
above follows since ni and 722 are arbitrary non-negative integers. The multiplicity of h{n + 1) is given by 
the number of pairs (?ii, 722) such that ni + n2 = n, which is precisely n + 1. 



'for instance Weyl quantization, but there are otlier possible choices 
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4.2 Quantization of the space S"^ x M^ and the Hamiltonians J and L 

We define the quantization of 3"^,^ to be the finite dimensional Hilbert space T-Le '■= ker(L — E). When 
E = h{n + 1), then dim(^£;) = n + 1 (otherwise T-Le = {0}). It will be convenient to introduce the 
"anihilation operators" Oj := -^(h-£r- + Xj), i = 1,2, which naturally quantize Zi/V^, i = 1,2 

respectively. Then L = h{aial + 0202 — 1). The quantization of the Hamiltonians x, y, z on 8^,2 are the 
restrictions to "H^; of the operators: 

X := -{aia*2 + a2al), y ■= — (aial - a2al) , z := -{aial - a2a*2). (4.3) 

This definition makes sense because TIe is stable under the action of x, y, z. This can be checked right 
away using the commutation relations [a^, a!|] = 1, but it will also follow from the explicit action of these 
operators, as explained in Section 1431 below. 

Of course, in M? ■., the quantization oivisv := (f^) and the quantization u of n is the multiplication 
by u (that we simply denote by u). Thus we have the very natural definition: 

Definition 4.1 The quantization of S'^,2 x M^ is the (infinite dimensional) Hilbert space T-Le ® L^(R) C 

L'^(M?) (g) L^(M). The quantization of J is the operator J = Id(g)[ - ^-^ + ^j + {z (g) Id). The 

quantization of H is the operator H = ^{x (Si u + y (g (j-^)). 

This definition depends on the energy E, which will be fixed throughout the paper. For the numerical 
computations, we have taken E = 2, which corresponds to the quantization of the standard sphere x^ + y^ + 

z^ = l. 

Lemma 4.2. The operators H and J commute, i.e. we have the identity [-ff , J] = 0, both in the functional 
analysis sense fie. as an unbounded operator on a dense domain), and in the algebraic sense, as a bracket 
in the Lie algebra of polynomial differential operators. 

Proof. It is enough to show that [iJ, J] = holds on elements of the form f ® g, where / is any element in 
Ue, and g G Cg°(M). And indeed, 

[H,J]{f®g) = {HJ-JH){f®g) = HJ{f®g)-JH{f®g) 
= H{f®Ng + {zf)®g)--{xf®ug + yf®vg) 

— 7T i^f ® uNg + xzf ®ug + yf® vNg + ij^f ® vg) 

— — {xf ® Nug + yf <S Nvg + zxf <S ug + zyf (8) vg) 

= xf(g)[u,N]g + [x,z]f(Sug + yf^[v,N]g + [y,z]f^vg. (4.4) 

As before, we have denoted N := —^-§;p: + \- Now 

l'*''"l^ = ''(-Td-^ + T)^-(-TS? + T)"^ = T(-«d-^ + d-^'')^ 

and 

^(uf) = /— + 2^— + n^ = 2^ + n^ 
du^ du^ dn dn dn^ du dv? 
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Hence [u, N]f = ^(2^) = h^^if)- Therefore [u, N] = ihv. Similarly, [v,N] = -ihu. It is also 
standard to check that the "angular momentum variables" {x, y, z) satisfy : [y, z] = —ihx and [x, z] = ihy. 
Hence expression (14.41 ) equals 

xf (g) {ihv)g + {ihy)f (S> ug + yf (S> {—ihu)g + {—ihx)f (E> vg = 0. 

The result follows. D 



Remark 4.3 Although the proof of Lemma 14.21 is interesting on its own, there is a theoretical reason for 
this lemma to be true, because our operators all derive from Weyl quantization of polynomial. And for 
such operators the following result is known: suppose that Hi is a quadratic Hamiltonian and H2 is any 
polynomial Hamiltonian function such that {Hi, H2} = 0. Then Moyal's formula |[26l ITTl [121 yields, 
formally, [i^i, H2\ = 0. In our case J is quadratic in the variables {u,v,xi,X2,Ci,C2)- This gives an 
alternative proof of Lemma |42l 



4.3 Joint spectrum of J, H 

We have left to find the spectrum of H and of J. First, we conjugate by the unitary transform in L^(M^) 

U : f{xi,X2) — )■ Vhf{Vhxi, vhx2)- 
This has the effect of setting h = lin the operator aj : 



]_(_d_ 
V2 \dxj 



U<^jU* = ^i— + xA=:Aj. 



Next, it is convenient to use the Bargmann representation 121, which states that the operator Aj defined 



a 



above and its adjoint A* are unitarily equivalent to the operators -^ and Zj, respectively, acting on the 

Hilbert space of holomoiphic functions on two variables L^p[(C^, 7r~^e~l^l ). (The notation Zj here is not 
exactly the same as the initial one in section l4~n but we keep it for simplicity.) 
The following lemma is standard. 

Lemma 4.4 (HI). The function ^^ ^ ^ , = -^, where a = (ai, 02), is an eigenfunction of L with norm 1 



and eigenvalue h{ai + a2 + 1). 

Proof. The function z°' is an eigenfunction of Zi-^ with eigenvalue Oj. Since L = h{zi-^ + Z2-^ + 1). 

we get L(z") = h{ai + 02 + 1)2;". 

We can compute ||z"|l , ,^„ , , ,2, = a!. Therefore the function 4^ is a normalized eigenfunction 

ofL. D 

Next we find the eigenspace of L for the eigenvalue h{n + 1). Since the monomials {z°^ / ^fa^-} a&i'^ 
form a Hilbert basis of the Bargmann space, the space T-Le = ker(I/ — h{n + 1)) is simply given by 

T-Le = span{^^ \ai+ a2 = n}, 
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thus it is the space of homogeneous polynomials of degree n in C^. We will use for it the following basis : 

\Z2 , Z1Z2 , . . . , Z-^ Z2, Zx J-. 

In order to understand the operator H, we need to consider z and N. The restriction of the operator 



k ^n—k\ 



/ -n 




2-n 








4-n 



z = f (aia^ — 020,2) to the Hilbert space He in given in terms of this polynomial basis by ^(^^1 2^2 

^(fc — (n — k))ziZ2~ ' . It follows that the matrix of f = ^{aial — 0202) relative to this basis is the diagonal 

matrix 

0\ 





h 
2 



V ... On/ 

Notice that this shows that He is indeed invariant under the action of z. Of course, a similar calculation can 
be done for x and y (see the proof of Proposition 14.7 I below) . Notice also that the eigenvalues of z range from 
— T^n to |n; in the case of the standard sphere S^ (with E = 2), we have the relation E = 2 = h{n + 1). 
Therefore the eigenvalues of z range from — ^ to -^. In the semiclassical limit n — )■ 00, we recover the 

=■ ° n+l n+1 ' 



classical range [—1,1] of the hamiltonian z on 5^. 

Next we consider the Bargmann representation for N 
L^„i(C^, TT-^e-l^l') and we obtain N = h{T-§p + i). 

The eigenf unctions of N are -^ corresponding to the eigenvalue h{k + ^] 

Lemma 4.5. The spectrum of J is discrete, and we have 

1 — 77, 



This time, we act of the Hilbert space 



spec (J) = h 



More precisely, for a fixed value X G ^(— 0^ + N), let £ 



2 



+ N 



-X := 



£x =spaniT^ ^ Z1Z2 ^ 



m + \ + k-l 



ker( J — A). Then 

A; < /c < n; 



>0 



In particular £\ has dimension 1 + min(r7, j + 



n-l^ 



Proof. In the double Bargmann representation, we have 

J 
Hence a simple computation gives 

J{t^ O zlz^-'') = h 



T 1 /^/ d 1,, h.d d . 

Id^(n(r- + -)) + -(z,--.2^) 



1 



H \-k 

2 



||(^' 



yk n- 
Z1Z2 



)Id. 



1 



(4.5) 



so the corresponding eigenvalues ai-e h{i + ^ + k — ^) where < A; < 77 and n, £ > 0. This shows that J 
admits a complete set of eigenvectors. Hence ker( J — A) is spanned by the set of eigenvectors coming from 
this family and corresponding to the eigenvalue A. This space is finite dimensional (hence J has discrete 
spectrum), and its dimension is the number of solutions {k,£) to the equation h{£ + | + fc — ^) = A with 
constraints < k < n; £ > 0, which is precisely 1 + min(r7, -^ + ^^). D 
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The fact that £\ is finite dimensional should be compared to the fact that the classical hamiltonian J is 
proper. 

Corollary 4.6. Given any n G N, and any A G ^(^^ + N), the ordered set 



r ^ 






A n 1 
fe = 0, 1, . . . ,mm(n, "^ + 2 ~ 2^' ^""^ 



Are 1 -1 



/£! v^A;!(re-A;)! 

/^ an orthonormal basis of £x- 

Our next goal is to compute the matrix of H. More precisely, since H commutes with J, the eigenspace 
£\ is stable by H. Thus, the spectral theory of H is merely reduced to the study of the restriction of H to 
£\, which we explicitly compute below. Then the best way to depict the spectra of J and H is to display 
the joint spectrum (see figure RTTT ). which is the set of (A, u) e M? such that, for a common eigenfunction /, 
one has both 

Jf = Xf and Hf = uf. 

Let £0 := f + f — |, ^ = mm{io,n) and let 




Figure 4.1: Semiclassical joint spectrum of J, H and momentum map image juxtaposed, computed using 
a numerical diagonalization of the band matrix in Proposition 14.71 In all our computations we have chosen 
E = 2, which corresponds to the quantization of the standard sphere x^ + y^ + z^ = 1. This implies the 
relation 2 = h{n + 1). Here re = 13, so /i ~ 1.14. 



Proposition 4.7. The matrix Mb;^ {H) of the self-adjoint operator H on the basis Bx is the symmetric matrix 



Mb AH) 





f ° 


Pi 






\ 




' /3i 





02 












0?. 





As 





h\i 
2) 










/3. 




I 









/3m J 
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Proof. We start by evaluating x and y on this basis: 

y{4^r') = j^ikzt'zr'^')-{n-k)z',+hr'^') 

We introduce: 

1 , ,dh ^ 1 dh. 

V2h^ du" y/2h^ du' 



Hence u{= u) = {a + a*)w |. Now we do the Bargmann representation 

h. d. ^ hd ia-Q*) fh 1 [h,d . 

2^^ +57^' " = T5^ = ^^V 2 = TV 2^37 - ^)- 



Hence we obtain 



n(/) = j;^(r^+i+^/-i), vi^r') = -J 'J^Ur'-' - r'+'-^ 



In what follows, for brevity of the notation, we write c^ := zJ^Zg"'^. Note that n is fixed. Recalling 

H = i(ii (g) X + u (8) y), we get 

- Qy^\^^'~' - r'^')ikck-i - (n - k)ck+i)) 

= \ (1)'^' [kr'^^ck-i + tt/-icfc_i + (n - A:)/+icfc+i + £(n - k)T'-^cu+i 

- ekT^~^Ck-i + i{n - ky~^Ck+i + kr^-^^ck-i - (n - ky^^Ck+i 



J)'^V/+icfc_i + (n - A;)^r^-icfe+i). (4.6) 

Notice how this formula, together with Lemma 1431 confirms that £x is stable under H. 

In order to have a better numerically prepared matrix (and a nicer-looking formula !), we next express 

everything in an orthonormal basis. Denote eik = -777 /^ ^ so that e^ k is an eigenvector of J of norm 

1: 



1 n 

J{ei^k) = h{i + - + k- -)ei^k = ^ei^k 






On the other hand we have that e^+i fc-i = , and that the first term of (14.71 ) is 



r^+^c._i 



V^£!A;!(n-A:)! ^£!A;!(n - fc)! 



^(£ + l)!(A;-l)!(n-A; + l)!e,+i,fc_i 



V(£ + l)fc(n-A: + l)ef+i,fc_i. 
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Similarly the second term of (14.71 ) is 

e{n - k)T^-^Ck+i i{n - k) 



yj{l-l)\{k + l)\{n-k-l)\ei_i. 



fc+i 



^l\k\{n-k)\ yj l\k\{n - k)\ 

= \Ji{k + l)(n - A;)e£_i,A:+i. 

Since ^ = ^o — ^. we get 
^(e£,fc) = (2) ^ (V(4 - fc + l)A:(n - A: - l)e^+i,fc_i + ^(^ - /c)(/c + l)(n - A:)e^_i,fc+i) 

//in, 3/2 

= (-^j (/3fce^+i,fc-i + /3fc+ie£-i,fc+i)- 
This, of course, gives the statement of the proposition. D 

4.4 The spectrum S(?t,) ^^ ii\^cx:U -\&) 

In the next section, we will be particularly interested in the J-eigenvalue A = 1, which corresponds to the 
J-critical value of the focus-focus point, in the classical system. Since E = 2 = h{n + 1), we see that 
(.Q = ^^^^ + ^^ = n. Therefore the dimension of ker( J — Id) is equal to ?i + 1. Notice that, for A < 1, 
the dimension of ker( J — A) is increasing linearly with slope 1 (with respect to the parameter k that we 
introduced above) whereas for A > 1 this dimension is constant, equal to ?i + 1. This can be seen as a 
quantum manifestation of the Duistermaat-Heckmann formula ||9l- 

5 Inverse spectral theory for quantum spin-oscillators 

The theme of this section is to give evidence of the following conjecture being true in the case of coupled 
spin oscillators: 

Conjecture 5.1. A semitoric system is determined up to symplectic equivalence by its semiclassical joint 
spectrum fi.e. the set of points in M^ where on the x-axis we have the eigenvalues A of J, and on the vertical 
axes the eigenvalues of H restricted to the \-eigenspace of J). From any such spectrum one can construct 
explicitly the associated semitoric system. 

In this section we try to convey some ideas to explicitly compute all the symplectic invariants from the 
semiclassical spectrum. It might not necessarily be the optimal way to prove an inverse spectral result, as 
some quantities are more easily defined implicitly rather than explicitly by the spectrum. But we believe 
that, from a quantum viewpoint, having constructive formulas for the symplectic invariants is particularly 
valuable. 

We emphasize the word "semiclassical" here : in order to recover the symplectic invariants we need be 
able to compute the joint spectrum for small values of h. What can be said for a unique, fixed value of h is 
much harder question. 

5.1 Polygon and height invariant 

Recovering the polygon invariant is probably the easiest and most pictorial procedure, as long as one stays on 
a heuristic level. Making the heuristic rigorous should be possible along the lines of the toric case explained 
in mi and lUll, but we don't attempt to do it here. 
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The first thing to do is to recover the image of the classical moment map, including the position of the 
singular values. This could be done by a local examination of density of the joint eigenvalues. 

Next, in order to recover the polygon invariant, we need to obtain the integral affine structure of the 
image of the momentum map. We know from ID |24l that the joint spectrum possesses a semiclassical 
integral affine structure on the regular values of the momentum map. This integral affine structure can be 
extended to the elliptic boundaries, as explained in |[24ll . Thus, except along a vertical cut through the focus- 
focus critical value, one can develop this affine structure such that the joint eigenvalues become elements of 
the lattice hZ^. See figure EU 




Figure 5.1: Recovering the polygon invariant. The top picture is the joint spectrum of {J,H). In the 
bottom picture, we have developed the joint eigenvalues into a regular lattice. One can easily check on this 
illustration that the number of eigenvalues in each vertical line in the same in both pictures. 

The convex hull of the resulting set is a rational, convex polygonal set, depending on h. Since the 
semiclassical affine structure is an /i-deformation of the classical affine structure, we see that, as ?i — ;■ 0, this 
polygonal set converges to the semitoric polygon invariant. 



5.2 Semiclassical formula for the spectrum E(n) 

In order to recover the Taylor series invariant from the spectrum, we need a precise description of this 
spectrum. There are two options : either describe the spectrum in regular regions, and then take the limit to 
the focus-focus critical value; or describe the spectrum directly in a small neighborhood of the focus-focus 
value. We choose the second option, because it seems more appropriate for a reasonably accurate numerical 
formula for the invariants, in the spirit of equation (12.281 ). 

The drawback of this approach is that there is no result currently available giving the description of this 
spectrum. The singular Bohr-Sommerfeld rules of |f22| would give the required result, in case J and H were 
pseudodifferential operators. Of course they aie not, since the phase space S"^ x M^ is not a cotangent bundle. 
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However they ai^e semiclassical Toeplitz operators, in the sense of H, and it is known that the algebra of 
Toeplitz operators is microlocally equivalent to the algebra of pseudodifferential operators 131. Therefore, 
we propose the following conjecture. 

Conjecture 5.2. The formula in Corollary 6.8 in Vu Ngoc's paper I122\I holds also if the operators therein 
involved are Toeplitz instead of pseudodifferential. 

This conjecture may be stated in the following way. Let S(n) be the spectrum of -f^lkcrfJ-idV ^'^^ 
bounded t G M, the formula 

A(t) - e(t)ln(2^) - 2argr(^^^^^^i^) G InZ + 0{ff^) 



holds if and only if ht £ S(n) + 0(^1°°) with 

(a) A(t) = X{t; h) admits an asymptotic expansion on integer > — 1 powers of h with smooth (=C°°) 
coefficients in t starting with \{t) = \ L "o + 170 {^(t)) + ^f + 0{h). 

(b) e(t) = e(t; h) has an asymptotic expansion on integer > powers of h with smooth coefficients in t 
starting with the second component of the vector .6(0, t) + 0{h) where B is the 2x2 matrix such 

thatB{J",H")^ = {qi,q2). 

(c) I^g{k{t)) is what is called the "principal value integral" of K{t), where i<,{t) is the 1-form on Aq 
defined by 

{k{t){Xj), k{t){XH)) = (0, t) ^^ {K{t){X,,), K{t)X,,)) = Bio, t) (5.1) 

Finally, /.^^ (k*) is defined in Proposition 6.15 of ll22l as 

0^Bo=7o{l-S2) s 

k{t) + e{t)ln{rAoPBo)) 
Ao=7o(si) ^ 

where e(t) is the first order term of e(t). 
For a semitoric system, the matrix B is of the form B = { ) , with B22 7^ 0. Thus we get 

V-D21 B22J 

e(t) = B22t. 
Moreover, because of formula (15.11 ). 

{nmXq,), R{t){X,,)) = (0, B22t). 

Therefore we see that -^^ = i?22K2,o. where K2,o is the restriction to Aq of the 1-form defined in equa- 
tion (12.31) . Thus, in view of equation (I2.4I ). we get an explicit formula for the symplectic invariant 02 : 

Though we haven't worked it out here, a similar- formula for the first invariant ai could be obtained 
along the same lines. 

In the case of the coupled spin-oscillator, B = i 9 ) > so B22 = 2 and 02 = ^^(^70 ('^*)) \t=o- 
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5.3 Obtaining 02 from the spectrum S(n) 

We show in this paragraph how the conjecture gives a way to obtain 02- Using formula (15.21 ) above, an easy 
corollary of the conjecture is Theorem 7.6 in [22], which says that 

. (Ek+i-Eu\ 27r/B22 , .^(^. ,, o. 

mm = 7- — ; h Oin) (5.3) 

V h y I In /il + a2 + In 2 + 7 ^ ' ^ ' 

for S(n) = {£^0 < -El < • • • < En}- Here 7 is Euler's constant. 

From the spectrum we can calculate t'^™{h) = min f ''+^~ ^ ) so 

B22{\\^h\+a2+\n2 + ^){l + 0{h)) = B22{\\^h\ + a2+\n2 + -^) + 0{h\nh). 



Therefore we may recover B22 as 



B22 = lim I . ^ ^, 1 . (5.4) 



Because the convergence of this Umit is very slow (of order |In h\ ), it is in practice much better to solve 
the system obtained with two different values of h, which gives : 



27r 2tt 



B22= *"""i%^^ /'"'"^'^^ +0{h^\nh^) + 0{h2\nh2). (5.5) 

Thus, if we choose h2 to be a fixed multiple of h = hi, we get a convergence speed of order 0{h\n h), 
which is indeed much more reasonable. 

Once B22 is known, it is easy to recover 02 > again through formula (15.31 ) : 

02 = lim I ^ - lln nl - In 2 - 7 ) , (5.6) 

and the convergence rate is again of order 0{hlnh). 

5.4 Numerical approximation of 02 using Maple 

Using Proposition |4.7[ we compute the spectrum S(n) of the Spin-Oscillator example for various values 
of n = 2/h — 1 by entering the matrix in the computer algebra system 'Maple' and ask for a numeric 
diagonalization. Then is it easy to implement the formulas ( 15.51 ) and (15.61 ). 

From the general theory, the minimal eigenvalue spacing is obtained — at least in the limit ?i — )• 0, at the 
focus-focus critical value H = 0. This is confirmed from the numerics. In fact, using the recursion formula 
for the characteristic polynomial Dn{X) of the matrix M^Bx {^) (with £q = n) : 

Dn{X) = XDn-l{X) - PlDn-2{X), 

we prove by induction that Dn{X) has the parity of n + 1. In particular, the spectrum is symmetric : 
S(n) = — S(n). When n is odd, is not an eigenvalue (D„(0) = (— l)("~^)/^/3i/33 • • • /3n), and hence the 
smallest spacing is simply twice the smallest positive eigenvalue : 

f^'^'in) = 2£[n]+2//i with h ^ 



n + l 
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Figure 5.2: Recovering the coefficient B22 (which is equal to 2 in our example). The horizontal scale is 
logarithmic: the integer abscissa k corresponds to n = 2^^ + 1. Thus h starts at 0.5 and decreases to the 
right to reach 1/513 ~ 0.002. The top curve — with circles — is the result of formula (15.41 ). which indeed 
converges very slowly. The curve with diamonds is obtained by the accelerated formula (15. 5I ). 
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Figure 5.3: Recovering the invariant 02. The graph plots the values of 02/ In 2 (which should be 5 in our 
example) computed using the formula (15.61 ). The horizontal scale is the same is in figure [ 



The results of our numerical experiments are plotted in figures |5^ and |5.3l They should be compared to 
the theoretical values of Theorem 12. 141 
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